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ABSTRACT 

We use abundances of Ca, O, Na, Al from high resolution U VES spectra of 2 00 red giants in 17 
globular clusters (GCs) to investigate the correlation found by iLee et al~l (|2009f ) between chemical 
enrichment from SN II and star-to-star variations in light elements in GC stars. We find that (i) the 
[Ca/H] variations between first and second generation stars are tiny in most GCs (~ 0.02 — 0.03 dex, 
comparable with typical observational errors). In addition, (ii) us i ng a large sample of red giants 
in M 4 with abundances from UVES spectra from iMarino et al. I (|2008[) , we find that Ca and Fe 
abundances in the two populations of Na-poor and Na-rich stars are identical. These facts suggest 
that the separati on seen in color-m agnitude diagrams using the U band or hk index (as observed 
in NGC 1851 by lHan et al. I l2009f ) are not due to Ca variations. Small differences in [Ca/H] as 
associated to hk variations might be due to a small systematic effect in abundance analysis, because 
most O-poor/Na-rich (He-rich) stars have slightly larger [Fe/H] (by 0.027 dex on average, due to 
decreased H in the ratio) than first generation stars and are then located at redder positions in the 
V, hk plane. While a few GCs (M 54, u> Cen, M 22, maybe even NGC 1851) do actually show various 
degree of metallicity spread, our findings eliminate the need of a close link between the enrichment 
by core-collapse SNe with the mechanism responsible for the Na-0 anticorrelation. 

Subject headings: stars: abundances — stars: evolution — globular clusters: general 



1. INTRODUCTION 

Stars of different stellar generations^ are currently rou- 
tinely found in all Galactic globular clusters (GCs). A 
common misunderstanding is to identify multiple se- 
quences in the color-magnitude diagrams (CMDs) as 
the only evidence of multiple stellar populations in GCs, 
while they are only the tip of the iceberg, observable 
when differences in chemical composition and/or age are 
very large. However, all GCs display large spreads in 
abundances of light elements. These arc due to the pres- 
ence of two generations of stars, separated by a small 
difference in age, not enough to appear as a smearing 
of the cluster turn-off, but clearly visible as a different 
chemical signatur e. This difference can be uncovered by 
spectroscopy (see iGratton et al. 1 120041 for a recent re- 
view). Since its dis covery by th e Lick- Texas group (see 
the early review by iKraft I [1994) , the Na-0 anticorrela- 
tion in GCs offers a powerful instrument to resolve age 
differences as small as a few 10 7 yrs (should very mas- 
sive stars be involved) by means of a signal as large as 1 
full dex in abundance. Whenever the Na-0 anticorrela- 
tion is observed in a GC, multiple stellar generations are 
present in the cluster. 

1 Based on data collected at the European Southern Observa- 
tory, Chile, programmes 072.D-507,073.D-0211,072.D-0742,077.D- 
0182 

2 INAF, Osscrvatorio Astronomico di Bologna, via Ranzani 
1, 40127, Bologna, Italy. eugcnio.carretta@oabo.inaf.it, an- 
gela.bragaglia@oabo.inaf.it, michelc.bcllazzini@oabo.inaf.it 

3 INAF, Osservatorio Astronomico di Padova, vi- 
colo dcll'Osscrvatorio 5, 35122 Padova, Italy. raf- 
faele. gratton@oapd . inaf . it ,sara.lucatello@oapd . inaf . it , 
valentina. dorazi@oapd .inaf. it 

4 Excellence Cluster Universe, Technischc Universitat 
Munchen,Boltzmannstr. 2, D-85748, Garching, Germany 

5 We will use the words stellar "generation" and "population" 
interchangeably, for the reasons explained later. 



The primordial P component of first generation stars, 
present in almost all Galactic GCs surveyed so far, 
has the high [O/Fe] and low [Na/Fe] ratios typical 
of core-collapse supernovae (SNe II) nucleosynthesis, 
characteristic of halo field st ars with simila r metal- 
licity. On the other hand, iCarretta et al. I (|2009al) 
showed that two other components of second gener- 
ation stars may be found: one dominant component 
has intermediate composition (I; observed in each GC), 
and another has extremely (E) modified composition, 
showing large O-depletion and Na (and Al) enhance- 
ments, present only in the most massive clusters^. 
This general pattern can be qualitatively produced by 
proton -capture reactions in H-burning a t high temper- 
ature (|Denisenkov fe Denisenkova 1 119891 : lLanger et al. I 
119931) . The fact that the same Na-0 and Mg-Al 
anticorrel ations are found in unevolved stars in sev- 
eral GCs (IGratton et al~ll2001l : lRamirez fe Cohen 1120021 : 
ICarretta et al. Il2004h tells us that this composition is in- 
herited from the ashes of a previous generation of stars. 

While the involved processes have been largely iden- 
tified, it is still not clear where these reactions oc- 
curred, the favorite sites being either main-sequence 
core H-burning of fas t rotating massive stars (FRMS, 
iDecressin et al. 1 12007( 1 or the hot bottom of the con- 
vective envelope in i ntermediate- mass asym ptotic giant 
branch (AGB) stars l| Ventura et al. 1120011) . Whichever 
the sites are, they can not synthesize iron or other el- 
ements heavier than Al. Indeed, apart from a few no- 
table except ions (uj Cen, see refe r ence in|Gratton et al. I 
[200l M 22. IMarino et al. ll200l t lDa Costa et al. Il2009fc 
M 54, ICarretta et al. I l2010aT Terzan 5. iFerraro et al. I 

6 The operational d istinction between I and E components 
llCarretta et al. Il2009al ) is not relevant here and we only consider 
the total fraction (I+E) of second generation stars. 
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[2001 and possibly NGC 1851. lYong fc Grundahl 112001 
the level in [Fe/ H] or a and Fe-group e lements is very 
constant in GCs (jCarretta et al. Il2009bl). 

However, a recent work by Lee et al*T()2009l . from now 
L09) seems to challenge this consolidated scenario, be- 
cause in their Ca-uvby survey they found a spread in the 
photometric hk index (including Ca n H and K lines) 
for several GCs. They interpret this spread as due to 
Ca abundance variations and claim that this a— element 
was produced by SN II in a past phase of the lifetime 
of GCs. In turn, these should be assumed to be initially 
much more massive than at present, else the highly ener- 
getic SNe winds would have been lost from the potential 
well. Moreover, they find that the Ca-strong group in 
GCs with evidence of multiple stellar populations from 
spectroscopy is associated to the population of stars with 
lower O and higher Na values, and vice versa concerning 
the Ca-weak stars. 

The need of massive proto-clusters is common to most 
scenarios of c hemical evolution of GC s, including one pro- 
posed by us (jCarretta et al. Il2010b|) . However a corre- 
lation between Ca and light elements would imply that 
the second generation stars were formed from material 
enriched not only by FRMS or AGB stars, but also by 
SN II events, which is a substantial modification of the 
current scenario. 

In the present Letter we exploit the large wealth of data 
from our ongoing proj ect "Na-0 anticorrclation and HB" 
(jCarretta et al. II2006I ) and from other extensive spectro- 
scopic surveys to compare [Ca/H] ratios obtained from 
high resolution spectroscopy with the abundance distri- 
bution of O, Na, Al in a sample of several hundred red 
giant branch (RGB) stars in 17 GCs. This allows to shed 
light on the role of Ca in the chemical evolution of GCs; 
from our extensive spectroscopic database, we seek to 
confirm whether the spread in the photometric hk index 
is driven by a real spread in Ca abundance. 

2. ATMOSPHERIC PARAMETERS AND ABUNDANCE 
ANALYSIS 

The full description of the analysis is given in 
ICarretta et al. I (j2009bl ): here we present only the abun- 
dances of Ca, which are based on a large number of 
lines (typically 10-12 in the most metal-poor clusters like 
M 15 and almost 20 in metal-rich GCs like 47 Tuc). The 
[Ca/Fe] ratios with the number of lines and the line-to- 
line r.m.s. scatter for each star are listed in Table [1] 
(completely available on line only) . In Table 2 we report 
the number of stars in each cluster, the average [Ca/Fe] 
and [Ca/H] values and the 1 a scatter of the mean. The 
rms scatters in [Ca/Fe] are smaller than in [Ca/H] be- 
cause the sensitivity of Ca and Fe lines to atmospheric 
parameters is quite similar. Five of the clusters listed in 
Table 1 are included in the sample by L09, who claim 
significant spread in Ca for four of them. 

The values listed in Table 2 indicate that the mean 
level of Ca is highly homogeneous in each GCs. The 
average r.m.s. star-to-star scatter for [Ca/Fe] is 0.03 dex 
from 17 GCs, to be compared with average scatter of 
0.18 dex, 0.22 dex, 0.07 dex and 0.25 dex for O, Na, 
Mg, and Al respectively. Note that the mean scatter in 
Mg, produced by core collapse SNe but also involved in 
proton-capture reactions, is more than twice that of Ca, 
which does not participate in any H-burning. 



3. RESULTS 

Since we have only about 10 stars per cluster observed 
with UVES, we combined the 17 GCs, separating first 
generation (P) stars from second generation (I+E) stars. 
For direct comparison with L09, we used [Ca/H] instead 
of [Ca/Fe], after normalizing the individual values to the 
mean of each cluster. In Figure 1 we show the histogram 
in [Ca/H] for P and I+E stars, with average values, er- 
rors, and rms. About two thirds of stars belong to the 
second generation, and the average normalized [Ca/H] 
values do not differ from those of first generation (at 
about la). A further check comes from the cumula- 
tive distributions, shown in Figure 2 (upper panel), from 
which we excluded the 12 stars of NGC 2808, the GC 
in o ur sample showing the most extreme var iations in 
He (jPiotto et al. I l2007t iBragaglia et~aTl 120101 and be- 
low in this Letter). The two distributions cannot be 
distinguished. Note that if we include also NGC 2808, 
the probability of the Kolmogorov-Smirnov drops to 
0.07, still compatible with no significant difference, but 
possibly indicating some actual variations in this par- 
ticular cluster. However, when we show in the three 
lower panels of Figure 2 the cumulative distributions for 
[O/Fe], [Na/Fe], and [Al/Fe], all elements involved in the 
(anti)correlations in GCs, we clearly see a striking differ- 
ence (compare e.g. to supplementary Fig. 14 of L09). 

Our sample of 200 stars comes from 17 different GCs, 
but there is already available in literatu re a sample of 
simila r size for a single GC, M 4, for which [Marino et al. I 
(|2008f ) analyzed high resolution UVES spectra of 105 
RGB stars. This cluster is also in the L09 sample; from 
their supplementary Fig. 14 (panel i) we see that almost 
all stars defined Ca-strong are Na-rich (in our interpre- 
tation, second generation I, E) and vice versa for the Ca- 
weak. However, when we divide the stars in the Marino 
et al. sample in P and I+E (using their same separation 
at [Na/Fe]=0.2) and compare the histograms and the 
means of the two populations (see Figure 3, left panel), 
they look very similaiO]. This is supported by the cu- 
mulative distributions (Figure 3, right panel): they are 
indistinguishable, according to the KS test. A further 
check was done considering only stars in the two extreme 
quartiles in Na abundance: the average Ca values in the 
more Na-poor and in the more Na-rich stars perfectly 
agree with the ones in Figure 3 for the I+E and P stars, 
respectively. 

The stars called Ca-weak (Ca-strong) on the basis of 
their hk index do not seem to have lower (higher) [Ca/H] 
values from direct analysis of high resolution spectra; this 
demonstrates that the differences in the hk index must be 
produced by something other than Ca variations. 

It is important to stress that the observational effect 
at the basis of the L09 claim corresponds to la spreads 
in [Ca/H] of about 0.03 dex in the most extreme cases 
(NGC 1851, NGC 2808, excluding the special cases of to 
Cen and M 22 where the presence of a significant spread 
in iron has already been established with high resolution 
spectroscopy) and of about 0.02 dex for typical cases (like 
M 4, M 5, and NGC 6752, see Table 20. In our view, 
these numbers call for two basic considerations: 

7 As noticed by the referee, the same result hol ds using the 38 
stars in NGC 6752 analyzed bv lYong et aT~l JUS). 

8 We used their supplementary Table 3 for the FWHM of the 
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1. Such tiny spreads are similar in amplitude to sev- 
eral uncertainties that are known to affect the 
abundance analysis (like e,g, variations in the He 
content); they are usually neglected since the over- 
all uncertainty on single measures is typically larger 
than this. We do not exclude that there might be 
dispersions at these very low levels, but presently 
they can not be separated from the intrinsic errors. 

2. even if we regard these very small variations as real, 
some basic algebra indicates that, assuming a mass 
of 2 x 10 5 M© and a metallicity of Z = 1CT 3 for the 
gas cloud that gave origin to the cluster, a single SN 
II event would be sufficient to produce the required 
amount of Ca (and Fe). 

4. THE CASE FOR NGC 1851 

If the tiny variations detected in most GCs are not 
measuring intrinsic scatter in Ca, what might be the 
cause for the variations in the hk index? L09 discussed 
several issues to show that the only way to affect the hk 
index is a change in Ca. Specifically, they made a strong 
case for NGC 1851, where the few "Ca-strong" stars 
(with abundances f r om hi gh resolution spectroscopy by 
lYong fe GrundahTI (|2008f ) are segregated to the red- 
dest positions in the V, hk CMD. In a recent paper, 
I Han et aT~l (|2009() presented a split of the RGB in this 
cluster in the U,U — I plane, where the redder sequence 
is populated by Ca-strong stars. 

We simply note the striking similarity between the 
C MD by Han et a l. wit h the U,U — B CMD for M 4 
by IMarino et al~l (|2008l ). In both cases, the RGB is 
clearly split into two distinct sequences. However, the 
sa me dichotomy of RGB sequences in M 4 is explained 
by IMarino et al. 1 (|2008[ ) as the N enhancement affecting 
the U — B colors (through the CN and NH molecular 
bands within the U filter), so that Na-rich (and likely 
N rich) stars are separated from Na-Poor (and N-poor) 
stars in the CMDs. Since we demonstrated in Figure 3 
that in M 4 the [Ca/H] ratio is identical for the Na-poor 
and Na-rich stars, the same conclusion might be valid for 
NGC 1851. 

A further p iece o f the puzzle comes from 
Carretta et al. I (l2009al) . who confirmed results by 
Yong et al. I 1)20081 ) in NGC 6752: Na-poor stars define 



a narrow N-poor sequence all confined to the bluest 
part of the RGB in NGC 6752, whereas the Na-rich 
(N-enhanced) stars are spread out to the red part of the 
RGB, reaching large N abundances. Now, regardless 
of the nature of the polluters, the ejected matter is 
alw ays enriched in the m ain outcome of H-burning, 
He. IBragaglia et al.l (|2010t in preparation) used a large 
sample (about 1400 RGB stars in 19 GCs) to show that 
stars with higher He content should have on average 
slightly larger [Fe/H] values the difference between the 
E and P stars being 0.027 dex in [Fe/H], even if there is 
no difference in the overall metal content Z. 

These findings suggest that the correlation observed 
by L09 between the (small) dispersion in [Ca/H] and 
the (large) spread in hk may be related to variations 

RGB of the GC we have in common, converted them back to scatter 
in hk assuming a Gaussian distribution, and further converted to 
scatter in [Ca/H] and [Fe/H] using the sensitivity of Ca on the 
variation in hk (suppl. Sect. 4) and the relations in Sect. 3. 



in He and light elements only, with no fresh production 
of Fe and Ca. The sizable effect on the Kolmogorov- 
Smirnov test of the inclusion of NGC 2808 in the sample 
(Sect. 3) lends further support to this scenario, as this 
cluster see ms to display especially large variations in He 
abundance (|Piotto et al. 1 120071: IBragaglia et al.ll2010D . 

Hence, while it is quite possible that in NGC 1851 
some small intrinsic scat ter in iron does actually exist 
(lYong fe Grundahl |[2008h . our results do not support the 
idea that core collapse SNe contributed to the enrichment 
of second generation stars in most GCs, and indicate that 
the production of the proton-capture elements occurred 
in sites not able to synthesize Ca and Fe-peak elements. 

5. DISCUSSION AND CONCLUSIONS 

Lee and coworkers devised a scenario for chemical en- 
richment in clusters with multiple populations. It in- 
cludes two-phases. First, core collapse SNe from first 
generation stars inject in the intra-cluster medium ma- 
terial enriched with Fe and a— elements; at least part of 
this material is not lost by the GCs owing to their larger 
mass at early epochs. Second generation stars form 
from this enriched material, incorporating the ejecta of 
intermediate-mass AGB stars with longer evolutionary 
times. The second generation stars are then produced 
from matter enriched both in heavy and light elements. 

However, in our opinion, this scenario has to face sev- 
eral problems: 

• Since there is clear evidence of multiple popula- 
tions from spectroscopy in all GCs, except per- 
haps a few small ones wit h only a han dful of 
stars analyzed (e.g. Pal 12, iCohen I l200lt Ter 7, 
ISbordone et al. |[2007l : both belonging to the Sagit- 
tarius dwarf galaxy), any proposed scenario must 
account for the formation of the generality of GCs. 

• While a single core collapse SN might be enough to 
produce the claimed scatter in [Ca/H], thousands 
(or even more) FRMS or massive AGB stars arc 
required to produce the observed pattern for ele- 
ments involved in proton-capture processes. The 
impact of these two mechanisms can not be the 



• Core collapse SNe should produce a (variable) 
enhancement in the mean abundances for other 
a— elements. No significant variation has been ob- 
served for Mg and Si (apart from those explained 
by the Mg-Al cycle; see Fig. 3 inlCarretta et al.l 
l2009bTl . or Ti (IGratton et al. II2004D . 

• In GCs like M 4, a split in hk of 45a p (hk) (measure- 
ment errors for the populations, L09, supplemen- 
tary Fig. 10) between two populations translates 
in virtually no difference in the [Ca/H] ratios ob- 
tained from a homogeneous large sample of stars 
with high resolution spectra. 

• If FRMS are the polluters, the whole model would 
fail, since these stars release the polluted matter 
in very short time, while still in MS. Hence, the 
timcscale is even shorter than that for production 
of core collapse SNe. 
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• The model by L09, adapted to explain the dou- 
ble RGB seen in NGC 1851, results in a number 
ratio of first-to second generation stars equal to 
3 (75%/25%). From spectroscopy, on the other 
hand, it was found that the bul k of stars in GCs be- 
long to the second gen eration (jCohen et al. 112001 : 
iCarretta et al. 1 1200 9V), with a number ratio typi- 
cally of about 0.5 (0.33/0.66). We stress that this 
feature is not a prerogative of mono-metallic GCs, 
the same value for the pr imordial fraction is als o 
found in u Cen and M 54 (jCarretta et al. Il2010ah . 
the two most massive clusters in the Galaxy, both 
showing actual (large) dispersion in [Fc/H]: this 
finding demonstrates that the relative fractions of 
first and second generation stars are not related to 
the enrichment of heavy elements. 

• The r— process element europium, whose main pro- 
duction site is expected to be core-collapse SNe, is 
fairl y constant within all t he studied globular clus- 
ters (|Gratton et al. 112004ft . 

• Finally, one of the strongest objection to the L09 
model comes from the run of Al. With its yield ex- 
tremely dependent on the neutron excess (metal- 
licity), any variations in the content of Fe (or 
heavy element produced in SN II, like Ca) should 
be accompanied by clear changes in the primor- 
dial Al abundance. As shown by ICarretta et al. I 
()2009bf ) this is not always seen, in particular in 
clusters like M 4 where a bimodal distribution is 
claimed in the hk distribution. The Al distribu- 
tion is far from being bimodal in M 4, on the con- 
trary [Al/Fe] is quite constant and does not vary 
between first and second generation stars. This 
is clearly seen also in other GCs like NGC 6171, 
NGC 288, and NGC 6838 from our own data 
(ICarretta et al. Il2009bf) as well in the larger sample 
by iMarino et al. I (|2008| ) . This constancy, due to 
the secondary (metallicity-dependent) production 
of Al in SNe strongly precludes that any signifi- 
cant difference in heavy metal-content is provided 
by SNe to stars in different stellar generations in 
GCs. 

We conclude that (i) the spreads in [Ca/H] are very 
small; (ii) different distributions in O, Na, Al are clearly 
seen between first and second generation stars in all GCs 
analyzed so far, on the contrary no statistically signifi- 
cant difference is found regarding Ca abundances from 
high resolution spectroscopy; (iii) the small difference in 
Ca observed in some cases might be well explained with 



systematic second-order effects due to the fact that sec- 
ond generation stars are He-enhanced, and appear more 
metal-rich and hotter when analyzed with high resolution 
spectroscopy. 

At the moment, we cannot offer a viable alternative 
explanation for the dispersion in the hk index observed 
by Lee and collaborators, although this issue certainly 
requires further studies. 

Partial funding come from PRIN MIUR 2007, PRIN 
INAF 2007, and the DFG cluster of excellence "Origin 
and Structure of the Universe" . 
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Table 1 

[Ca/Fe] and [Ca/H] ratios for individual stars. Complete table available on line only. 



GC ID 


star 


nr. 


[Ca/Fe] 


a 


[Ca/H] 


flag PIE 


NGG0104 


5270 


15 


0.299 


0.150 


-0.473 


2 


NGG0104 


12272 


11 


0.310 


0.107 


-0.437 


2 


NGG0104 


13795 


11 


0.345 


0.105 


-0.487 


1 


NGG0104 


14583 


11 


0.309 


0.110 


-0.375 


2 


NGG0104 


17657 


11 


0.298 


0.127 


-0.546 


3 


NGG0104 


18623 


19 


0.322 


0.213 


-0.451 


2 



a nr is the number of Ca lines in each star. 

b flag is 1,2, and 3 for the primordial P, intermediate I and extreme E component, respectively; flag means that either Na or O was missing 
and the star could be no assigned to a component. 
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Table 2 

Mean [Ca/Fe], [Ca/H] ratios and internal errors. 



GC ID 


lu- 


[Ca/Fe] 


a 


errl 


err2 


[Ca/H] 


(7 


NGC 0104 


ll 


+0.31 


0.01 


0.03 


0.03 


-0.45 


0.06 


NGC 0288 


10 


+0.41 


0.02 


0.02 


0.03 


-0.90 


0.05 


NGC 1904 


10 


+0.28 


0.01 


0.02 


0.02 


-1.30 


0.04 


NGC 2808 


12 


+0.34 


0.02 


0.03 


0.03 


-0.81 


0.07 


NGC 3201 


13 


+0.30 


0.03 


0.02 


0.02 


-1.21 


0.07 


NGC 4590 


13 


+0.26 


0.03 


0.03 


0.03 


-2.00 


0.04 


MPP Z.QCIA 
ly\jKj OiJ\J"± 


14 


+0.38 


0.02 


0.02 


0.03 


—0.96 


0.05 


NGC 6121 


14 


+0.42 


0.03 


0.02 


0.03 


-0.75 


0.07 
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a nr is the number of stars in each cluster. 

errl for [Ca/Fe] account for internal errors in T B ff, 


Vt and EW , whereas err2 includes all 


sources of uncertainties 


(see Carretta et al. 


2009b). 
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Figure 1. Distribution of [Ca/H] ratios from high resolution UVES spectra for 202 RGB stars in 17 globular clusters. The cross-hatched 
red area is the histogram for stars of the primordial P component, while the blue hatched area indicates stars of the second generation 
(I+E stellar component, see text.) 
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Figure 2. Cumulative distributions for stars of the P component (solid red lines) and for second generation stars (I+E components, 
dotted blue lines). From top to bottom the distribution of abundances of [Ca/H], [O/Fc], [Na/Fc] and [Al/Fc] arc plotted. In the top panel 
NGC 2808 is excluded (see text); if included the KS probability would become 0.07 . 
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Figure 3. Distribution of [Ca/H] ratios in M 4 (from UVES spectra bv lMarino et a l, 2008; left panel), and cumulative distributions for 
stars more Na-rich (blue) and more Na-poor (red) than [Na/Fe]=+0.2 dex, that we identify with the I+E and P components, respectively. 



